Dynamical holographic QCD with area-law confinement and linear Regge trajectories 
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We construct a new solution of five-dimensional gravity coupled to a dilaton which encodes es- 
sential features of holographic QCD backgrounds dynamically. In particular, it implements linear 
confinement, i.e. the area law behavior of the Wilson loop, by means of a dynamically deformed 
anti-de Sitter metric. The predicted square masses of the light-flavored natural-parity mesons and 
their excitations lie on linear trajectories of approximately universal slope with respect to both 
radial and spin quantum numbers and are in satisfactory agreement with experimental data. 
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Over the past decade a qualitatively new perspective 
on strong- interaction physics emerged from gauge/string 
dualities 1] and the underlying holographic principle. 
These dualities map (i) string theory spectra on asymp- 
totically AdS xX spacetimes (i.e. AAdS xAT, where X 
is a compact space) into gauge invariant, local operators 
of the dual field theory, (ii) the fields parameterizing the 
boundary conditions into sources for the dual operators, 
and (iii) the string theory partition function (or its low- 
energy gravity limit) into the generating functional of 
the field-theory correlators. As a consequence, the noto- 
riously complex strong-coupling regime of large- iVc gauge 
theories can be approximated (in low-curvature regions) 
by weakly coupled and hence analytically treatable classi- 
cal gravities. The gauge/gravity correspondence thereby 
supplies new analytical tools for the study of hadronic 
observables in the non-perturbative regime of the strong 
force. 

Applications of gauge/gravity dualities to "QCD-like" 
gauge theories either start from specific D-brane setups 
in ten- (or five-) dimensional supergravity and derive the 
corresponding gauge theory properties, or try to guess a 
suitable background and to improve it in bottom-up fash- 
ion by comparing the predictions to QCD data. Even the 
simplest and oldest bottom- up (or "AdS/QCD") model, 
the hard wall reproduces a surprising amount of 
hadron phenomenology The conformal invariance of 
AdSs in the UV reproduces, in particular, the counting 
rules which govern the scaling behavior of hard QCD 
scattering amplitudes, while an infrared cutoff on the 
fifth dimension at the QCD scale Aqcd implements the 
mass gap and discrete hadron spectra. 

The hard-wall predictions for the squared masses of 
light-flavor hadrons depend quadratically on the princi- 
pal and spin excitation quantum numbers however, 
in contrast to the theoretically expected and approxi- 
mately observed linear Regge behavior A straight- 
forward way to correct this shortcoming was suggested 
m Ref. '^I where the AdSs geometry is kept intact while 



an additional dilaton background field with quadratic de- 
pendence on the extra dimension is exclusively respon- 
sible for conformal symmetry breaking. This dilaton 
soft- wall model indeed generates linear Regge trajectories 
TOj^ g ^ n-\- S for light-flavor mesons of spin S and radial 
excitation level n. (Regge behavior can alternatively be 
encoded via IR deformations of the AdSs metric @, I3]-) 

However, the dilaton soft wall fails a quintessential test 
for conflning gauge theories: the resulting vacuum expec- 
tation value (vev) of the Wilson loop does not exhibit the 
area-law behavior which a linearly conflning static quark- 
antiquark potential would generate. This is because the 
Wilson loop vev is determined by the area of the dual 
string world sheet in the five-dimensional spacetime 0], 
i.e. it depends exclusively on the background geometry. 
Since the latter remains exact AdSs (and thus conformal) 
in the soft-wall model, the Wilson loop vev shows a non- 
confining perimeter law. (The hard wall, on the other 
hand, also confines magnetic charges instead of screening 
them [5|.) A second, common shortcoming of both hard- 
and soft-wall backgrounds is that they are not solutions 
of a dual gravity. Hence their relation to the dynamics 
of the sought after QCD dual remains obscure and all 
gauge theory vacuum properties (including confinement, 
chiral symmetry breaking and condensates) have to be 
imposed by hand [13] ■ 

In the following we show how both shortcomings can 
be overcome, by deriving a rather minimal AdS/QCD 
background which implements the area law, i.e. lin- 
ear confinement, dynamically. Our strategy will be to 
adopt an IR-deformed (and hence non-conformal) AdSs 
ansatz family for the metric which is general enough to 
generate the area law while keeping the fifth dimension 
non-compact in order to allow for linear Regge trajecto- 
ries. We will then construct the corresponding dilaton 
fields and potentials such that their combination solves 
the five-dimensional Einstein-dilaton equations, and find 
those solutions which additionally generate linear Regge 
trajectories in the highly excited meson spectrum. 
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To implement the above program, we start from the 
Einstein-Hilbert action of five-dimensional gravity cou- 
pled to a dilaton <i>, 
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(1) 

where k is the five-dimensional Newton constant and V 
is a still general potential for the scalar field. We then 
search for static solutions of the corresponding field equa- 
tions in which the metric is restricted to the form 



gMN 



(2) 



with J] = diag(l, — 1, — 1, — 1, — 1) and the dilaton field 
$ (z) depends on the radial coordinate only. We write 
the warp factor as 



A{z) ^Inz + C'iz) 



(3) 



(our units are such that the AdSs curvature radius is 
unity) where the function C (z) describes non-conformal 
deformations of the AdSs metric. We further impose the 
boundary condition C(0) = which restricts the geome- 
try to asymptotically AdSs (AAdSs) spacetimes and thus 
ensures conformality in the ultraviolet (UV). 

Variation of the action ^ leads to the Einstein-dilaton 
equations 
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0, (4) 
0, (5) 

(6) 



for the background fields A and $. This set of coupled 
differential equations is redundant, i.e. only two of them 
are independent 0. To cast those into the form most 
suitable for our purposes, we add the two Einstein equa- 
tions to obtain the derivative of the dilaton field as a 
function of the warp factor. 



$' - VsVA'^ + A" 



(7) 



(the positive sign is chosen for definiteness). Equation ^ 
determines (up to a constant to be fixed by a boundary 
condition) the dilaton field which forms, in combination 
with a given warp factor A, a solution of Eqs. ([1]) - (O. 
Equation ([7]) shows, in particular, that a constant dilaton 
solution requires the warp factor to satisfy A" = —A''^, 
i.e. under C {0) — exact AdSs is the only solution. 
In other words, the conformal-symmetry breaking solu- 
tions of interest require a non-constant dilaton and conse- 
quently an IR deformation of the AAdSs geometry. (The 
lack of the latter in the soft-wall background ^ explains 
why it cannot solve Eqs. ([4]) - ([6]), as realized in Ref. 
[ll|.) Another consequence of Eq. ^ is that there re- 
mains only one boundary condition to be imposed on the 
dilaton field. 



In order to recover the full information content of the 
field equations (|4]) - ([6l), it remains to find the dilaton 
potential V (z)) (whose specific form played no role 
in the above discussion) for which a given A and the 
corresponding $ determined by Eq. ([7|) form a solution. 
This potential can be obtained either from one of the 
Einstein equations ([4]), ([5]) or (up to a constant) from 
the dilaton equation by substituting Eq. In 
either case, the result is the relation 



[A"(z)-3A'2(z)] 



(8) 



between the dilaton potential, evaluated at the dilaton 
solution, and the metric. With the help of Eq. ([8]) and 
^'dV/d^ — dV/dz it can immediately be verified that 
all field equations ([3]) - (O are indeed fulfilled by dilaton 
fields which satisfy Eq. Hence Eqs. ^ and © 

are equivalent to the two independent Einstein-dilaton 
equations, and they allow to obtain explicit solutions <I> 
and the corresponding potential {z)) from a given 
warp factor A (without requiring a superpotential) . 

Having shown how to construct such solutions, we can 
now pick a specific warp factor suitable for our purposes. 
To this end, we recall that the vacuum expectation value 
of the Wilson loop in the dual boundary theory is re- 
lated to the area of the string world sheet in the five- 
dimensional space-time 3] and thus determined by the 
geometry ([2]) alone. Hence we adopt a subfamily of met- 
rics which encodes (besides the conformal UV behavior) 
the desired area law behavior. The specific conditions 
given in Ref. have been used in Ref. to show 
that the leading IR (i.e. z oo) behavior of a linearly 
confining (and magnetically screening) metric of the form 
([2]) in a non-compact fifth dimension is characterized by 
C(z) = z'^ -I- . . . where A > ensures that the conformal 
AdSs metric dominates the UV limit (i.e. C (0) — 0) and 
A > 1 is required for linear confinement. 

In order to transparently elaborate on the impact of 
such confining C{z) on the meson spectrum, we first ne- 
glect subleading terms and begin our quantitative discus- 
sion with the minimal choice 



(9) 



for the confining AAdSs warp factor. We then obtain the 
corresponding dilaton solution by integrating Eq. ^ un- 
der the UV-conformality preserving boundary condition 
$(0) = 0. The result is 

$a(z) = ^ [(1 + A)ln (Az* + VA + A2 + A2za) 
- (1 + A) In (VA + A2) + z^/ VA + A2 + A2z^] (10) 

which, together with the AAdSs metric specified in Eq. 
(121), forms a new analytical solution family of the five- 
dimensional Einstein-dilaton equations (|31) - ([HI), with 
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the dilaton potential given by Eq. JS]). For A = 2 it 
reduces to a solution given in Ref. [9|. (To elevate the 
soft-wall dilaton field and the AdSs metric to a solution 
of the coupled field equations, in contrast, requires an 
additional tachyon field as suggested in Ref. [Hj.) 
Evaluated at the solution the dilaton potential is 



(11) 



Close to the UV and IR limits the dilaton solution be- 



haves as ^\{z) 
Co 



Cqz^/'^ and (^\{z) 
V3. 



CooZ" where 
In these limits 



2v/3(l + A-i) and Coo 
the potential can be easily expressed as a function of 
i.e. F($) -6 + 3(A + 1)(A - 8)$V (2c2) and 

F($) -9A2$2exp(2$/coo) / (2c^). The exponen- 

tial divergence for $ — > oo is both induced and coun- 
terbalanced by the exponential behavior of the metric. 
Indeed, the dilaton-gravity coupling is necessary to sta- 
bilize the classical solution since the dilaton potential is 
not bounded from below. As expected, the boundary 
condition <&(0) = ensures conformal invariance in the 
UV where the potential V^(0) = —6 consists of just the 
AdSs cosmological constant term. 

We will now derive the mesonic excitation spectrum 
in the background ([9]), (fTO)l and then study the existence 
criteria for another classic confinement signature, namely 
linear square mass trajectories of Regge type for highly 
excited hadrons. To this end, we utilize the tensor gauge- 
field framework and notation of Ref. Q . 

The spin S string modes of the massive tensor fields 
4>Mi...Ms (ill axial gauge) in the dilaton-gravity back- 
ground can be rewritten in terms of reduced amplitudes 
V'n.s which satisfy the Sturm-Liouville equation 



where the spin-dependent string-mode potential is 



Vs{z) = 



B'^{z) B"{z) 



(12) 



(13) 



with B ~ {2S — 1) A + ^. The gauge/gravity dictionary 
identifies the eigenvalues g with the squared meson 
mass spectrum of the boundary gauge theory. Equa- 
tion (fT5)) implies that the mode potential Vs depends 
on A and only, i.e. not on $ and its boundary con- 
dition. Since both the dilaton solution and its potential 
are determined by the confining ansatz for the warp fac- 
tor A (cf. Eqs. dZl) and ([8])) and by UV conformality, 
the mode potential may also be expressed exclusively in 
terms of the warp factor. Since our main focus is on the 
higher-spin mesons and their mass trajectories, further- 
more, we have derived the dual meson dynamics (|12p 
from the quadratic part of the action only, as in Ref. 
0. We have also neglected additional nonlinearities due 
to higher-dimensional operators which would arise from 



stringy a' corrections in the gravity action and in the 
non-Abelian DBI action of flavor brane stacks, or from 
interactions with additional bulk fields. Their neglect 
is common practice in the current, first generation of 
AdS/QCD duals. A more detailed discussion of the en- 
suing hmitations can be found e.g. in Refs. (5l.[loj. 

Important qualitative aspects of the meson spectrum 
arising from Eq. (I12p can be understood by studying the 
UV (i.e. z ^ 0) and IR {z oo) limits of the mode 
potential. We start with the UV behavior. Keeping only 
the leading terms for small z and A > 1, 

Vs{z) = ao (S) z-^ + ai {S) z^-^ + a2 {S) z^-^ + . . . (14) 

with the spin-dependent coefficients aQ{S) — 5*^ — 1/4 
(which originates from the AdSs part of the warp fac- 
tor), ai(S') = ^3A(A-f 1)(S'- A/4) and 02(6') = 
A [SS-^ - 4(A + 1)5" 5A 3] /4. For A = 4 the dilaton 
field therefore becomes proportional to z"^ when z — > 
(see the discussion below Eq. pTjl l. as imposed by hand 
in Ref. However, while in this case the 02 term in Eq. 
(|14p becomes harmonic and generates approximately lin- 
ear trajectories for the low-ly'mg radial excitations, the 
mode potential will grow as z^ at large z and cause strong 
nonlinearities in the trajectories of high-lying radial and 
spin excitations. 

To implement the linear trajectories for highly excited 
meson states, which semiclassical arguments predict as a 
consequence of long, unbroken flux tubes in the large- A^c 
limit, we now turn to the large-z behavior of the mode 
potential. The leading infrared contribution is 



V5 {z) 



A2 



{2S + Vi-lf z 



2A-2 



(15) 



Hence for A > 1 the entire normalizable spectrum is dis- 
crete and has a mass gap (in the absence of normalizable 
zero modes), as previously found for radial glueball exci- 
tations Q and in remarkable agreement with the area-law 
condition A > 1 Q for linear confinement. (In contrast, 
the soft-wall model generates a non-confining perimeter 

law :8g.) 

Asymptotically equidistant g values corresponding 
to linear trajectories further require A = 2 in the IR dom- 
inant part (l9|) of the warp factor, in order to generate a 
harmonic IR potential. The trajectories resulting from 
the simplest choice C{z) — z^ would be of non- universal 
slope (rn^ g oc n(S'-l-c)), however, which can be cor- 
rected by adopting 



C{z) = 



1 + a/3 (zAqcd)' 
25* + V3- 1 l + e(i-^AQCD) 



(16) 



as the non-conformal warp factor. The corresponding 
metric remains close to AdSs in the UV but deforms 
rather rapidly for z > Aqq^ to approach the confining 
large-2 asymptotics of Eq. © with A 2. The asso- 
ciated dilaton field and potential, which turn Eqs. ([31), 
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(116)) into a solution of the Einstein-dilaton equations (|4|) 
- ([6]), are then obtained by solving Eq. ([7]) numerically 
according to the procedure outlined above. This deter- 
mines the potential (|13p . and the masses follow by solving 
the eigenvalue problem numerically. In Fig. ^ the 
resulting spectrum is compared to experimental data and 
hard- and soft-wall model results. Note that the state la- 
bel n > 1 is chosen such that n = 1 refers to the nodeless 
ground state of the radial excitation spectrum. A sat- 
isfactory description of the meson mass spectrum with 
nearly universal Regge slopes is indeed achieved without 
any tuning of adjustable parameters. (The spin depen- 
dent factor in Eq. (|16p is required by universality. For a 
physical interpretation see Ref. A good analytical 

approximation to the spectrum for Aqcd — 0.3 GeV is 
(in units of GeV) 

m^^c. 1(11,^ + 95-9), in>l) (17) 



which makes the approximate universality of the linear 
trajectory slopes explicit. 




1 2 3 4 5 1 2 3 4 5 6 

n S 
FIG. 1: (a) Radial excitations of the rho meson in the hard- 
wall (dashed line), soft-wall (dotted line) and our dynam- 
ical soft-wall (solid line, for Aqcd = 0.3 GeV) backgrounds, 
(b) Square mass predictions of spin excitations compared 
to PDG 0] values for p(770), lj(782), $(1020), 7ri(1400), 
i.j(1420), p(1450), p(1570), 7ri(1600), 1.^(1650), $(1680) 
p(1700), p(1900), p(2150), /2(1270), a2(1320), /2(1430), 
/^(1525), /2(1565), /2(1640), a2(1700), /2(1810), /2(1910), 
/2(1950), /2(2010), /2(2150), /2(2300), /2(2340), c^3(1670), 
P3(1690), $3(1850), P3(1990), P3(2250), 04(2040), /4(2050), 
/4(2300), P5(2350), 06(2450) and /6(2510). 

Since the AdS/CFT dictionary translates the z depen- 
dence of the dilaton into the running of the gauge cou- 
pling, the latter could be implemented into Eq. (jl6[) 
for small z <C Aq^^ according to the perturbative 
QCD (3 function. This generates the leading correction 
Cpcrtiz) = -(21nz)-i [9, 11] which modifies the UV be- 
havior of the string mode potential (fT4)) as Vs {z) 

2-^ — \- \ t; n? , s -f . . . but vanishes in the IR. Hence 

asymptotic freedom and the perturbative corrections to 
it can naturally coexist with confinement at large z. 

To summarize, we have found a new solution of the 
five-dimensional Einstein-dilaton equations which pro- 
vides an approximate dual background for holographic 



QCD. The method used in its derivation applies to essen- 
tially all AAdSs (and hence UV conformal) spacetimes 
with a Poincare-invariant boundary. The vacuum proper- 
ties of the boundary gauge theory, including quark con- 
finement and condensates, are dynamically encoded in 
this solution without the need for additional background 
fields. In particular, our background generates a con- 
fining area law for the Wilson loop (in contrast to the 
soft- wall model), and we have outlined how the pertur- 
bative running of the gauge coupling could additionally 
be implemented. 

For metrics whose warp factors approach the power 
law z^ in the infrared of the non-compact extra dimen- 
sion, we found the existence condition for the confining 
area law, A > 1, to essentially coincide with the condition 
A > 1 for an entirely discrete meson spectrum and for the 
existence of a mass gap, as previously encountered in the 
glueball sector. With A = 2 our background solution sat- 
isfies this condition, and it organizes the square masses 
of both radial and spin excitations into linear trajecto- 
ries. It further generates the approximately universal 
slope of the experimentally observed trajectories (with- 
out introducing adjustable parameters beyond the QCD 
scale) and satisfactorily reproduces the empirical spec- 
trum of the light-flavored natural parity mesons. 
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